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Abstract 

An algorithm to measure the drift time of charge carriers in p-type point 
contact (PPC) high-purity germanium (HPGe) detectors from the signals 
processed with a charge-sensitive preamplifier is introduced. It is demon- 
strated that the drift times can be used to estimate the distance of charge 
depositions from the point contact and to characterize losses due to charge 
trapping. A correction for charge trapping effects over a wide range of ener- 
gies is implemented using the measured drift times and is shown to improve 
the energy resolution by up to 30%. 

Keywords: germanium detectors, neutrinoless double-beta decay, charge 
trapping 



1. Introduction 

P-type point contact (PPC) high-purity germanium (HPGe) detectors 
[H, lil have excellent properties for use in research in astro-particle and nuclear 
physics. They have generated substantial interests in neutrinoless double- 
beta decay searches by the Majorana Demonstrator Q and GERDA 
collaborations j3l due to their ability to distinguish between multiple and 
single-site interactions jsi [oj. The weighting potential in these detectors is 
sharply peaked near the point contact, giving rise to distinguishable current 
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pulses for simultaneous charge depositions separated by as little as 1 mm. 
Background from multiple-site charge depositions, such as gamma rays scat- 
tering in the detector, are thus distinguishable from neutrinoless double-beta 
decays which deposit their energy in a small volume. 

Another feature of these detectors is that the low capacitance of the small, 
point-like, collecting contact results in very low levels of electronic noise. This 
allows the use of PPC detectors with energy thresholds below 0.5 keV and 
has made them attractive candidates for direct dark matter detection [2I, 0]. 

The weaker fields in a PPC detector present some challenges. For ex- 
ample, determining the exact start time of the drift of charge carriers can 
be difficult as the signal generally originates in a region of weak weighting 
potential where it may be obscured by noise. The drift time is dependent 
on the location of the ionizing interaction in the crystal; measuring it can 
thus be a way to monitor the distribution of energy depositions, which could 
be of use in identifying sources of background in an experiment. Another 
challenge is the longer drift times in PPC detectors which can result in a 
degradation of the energy resolution due to charge trapping by impurities in 
the crystal. It has been shown that the drift time of events can be used 
to correct the measured energy, resulting in a substantial improvement in 
energy resolution. 

This paper introduces a digital algorithm to improve the determination 
of the drift time of events using the rising time of the charge pulses. The 
performance of this method is then investigated using simulated pulse shapes 
from energy depositions at different points in a detector. It is shown how 
the drift times can be used to obtain an estimate of the distance between 
charge depositions and the point contact and how this is approximated by 
their z-position (along the longitudinal axis of the crystal). Additionally, it 
is shown how the charge trapping in a detector can be examined and how the 
correction introduced in [sj can be improved by using an optimized measure 
of drift time. 

2. Pulse drift time measurements and validation 

PPC detectors have a sharply localized weighting potential around the 
point contact. Charges thus induce most of the signal on the electrode when 
they are closest to the point contact. An induced pulse will generally have 2 
components, a slow component when the charges are far from the electrode 
and a fast component when the charges are close to the electrode. Traditional 
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rise time metrics (such as the time between the 10% and 90% points of the 
charge pulse, as used in jsj) are thus biased towards measuring the time that 
charge carriers spend close to the electrode rather than the actual drift time. 

Figure [1] shows a comparison of the start time, tio, defined as the 10% 
rise of the charge pulse (solid) and, io.i, defined as the 0.1% rise of the 
average current pulse (dashed). Events of two different energies (691 keV 
and 58keV) are shown. For the higher energy event (panel (a)), the charge 
waveform could be used to determine an earlier start time; however, this 
would be much more difficult to do for the lower energy waveform because of 
noise. One should also note that, in both cases, using tio misses much of the 
slow component of the pulses, which is most directly correlated to the time 
that charges drift in the crystal. 




Figure 1: Charge pulse (solid line) with average current pulse (dashed line) for a 691 keV 
(panel a) and 58kcV (panel b) event. The vertical lines show, tio, the 10% rise time of 
the waveform (solid line) and, io.i, the 0.1% start time of the current (dashed line). The 
traditional tio misses much of the actual drift time in the initial slow component. Using 
the charge pulse to determine a start time is difficult because of noise, as seen in the lower 
energy event. The high energy event is a characteristic 'multi-site' event, likely from a 
gamma ray interacting twice in the crystal. 



The difficulty in determining an accurate pulse start time lies in removing 
biases introduced by noise. Different methods were explored for removing 
the noise without biasing the start time and it was found that the average 
derivative of the pulse resulted in the least bias (see Figure El discussed 
below). Methods to de-noise the waveform, such as a Daubechies wavelet- 
based filter or a running average were found to introduce distortions in the 
initial rising edge of the waveform that resulted in a small bias towards earlier 
start times. 
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A filtering algorithm that performs a running linear regression to a dis- 
cretely sampled signal, developed for various PPC detector pulse-processing 
applications, was found to be useful here in determining the average deriva- 
tive of a pulse for use in the drift time measurement. Given a set of samples 
from a digitized charge waveform, g^, collected from the output of a charge- 
sensitive preamplifier at discrete times ti and an averaging window length, 
W, the algorithm performs the linear regression, Oj + Sji, to the signal in the 
time window from i to i + W. This effectively gives the slope, Si, and the 
offset, Oi, of a linear fit to the signal between i and i + W. The slope and 
offset for each value of i were calculated as: 

i+W i+W 

k=i k=i 
i+W i+W 

k=i k=i 



W{xyi) - {x,){yi) 

(1) 



W{x])-{x,){x,) 
(Vi) - Si{xi) 



W 

The pulse defined as the set of Sj is the desired average derivative of the 
charge signal over W samples, and is defined herein as the average current 
pulse. Additionally, the algorithm also gives the average waveform, defined 
by the set of (yi), integrated over the same number of samples. Iq i, the point 
at which the averaged current reaches 0.1% of its maximum value (determined 
by moving backwards from the maximum) was found to be a good point for 
defining the start of the drift time. An optimal value for the averaging 
window was determined to be 500 ns in the apparatus discussed later in this 
paper. This was found to be most sensitive to the slow component of the 
charge pulse without introducing any significant bias. The end time of the 
drift, tgo, can be taken as the point where the charge pulse reaches 90% of 
its maximum. 

Figure [2] shows a comparison of the distributions of start times determined 
using the average current and the average waveform. In both cases, the 
start time was determined as the point when the average current (or average 
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waveform) decreased to 0.1% of its maximum value (stepping backwards 
in time), with a 500ns averaging window. The data collected correspond to 
events between 58 keV and 62 keV, where noise makes it difficult to determine 
the starting time of the drift (see Figure [1]). There is a clear bias for the 
average waveform-based calculation to give earlier rise times that did not 
correlate with the known position of the events in the crystal. 
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Figure 2: Comparison of the distributions of starting times determined as the 0.1% rise 
of the average current (solid) and average waveform (dashed). The average current and 
waveform were both determined using a 500 ns averaging window and the algorithm in 
equation [TJ The data are energy depositions between 58keV and 62kcV collected with 
the custom point contact detector described in section |4l There is an apparent bias for 
the waveform-based calculation to determine an (unphysical) earlier start time. 



3. Monte- Carlo validation of drift time measurement 

A point contact detector with the same geometry and charge impurity 
profile as the detector used in the measurements (section |4]) was simulated. 
The detector has a radius of 30.75 mm and a height of 50 mm. The simu- 
lation was based on code developed by the GRETINA collaboration joj and 



is described in [10|, lllj • Charge pulses from the detector were simulated for 
energy depositions on a 2 mm x 2 mm grid without including the electronics 
chain or noise. Gaussian noise, based on the amplitude of the pulse, was later 
added to the simulated pulses with amplitudes consistent with that seen in 
the data (1% noise for a 60keV event, 0.1% noise for a 1 MeV event, relative 
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to the maximum amplitude of the pulses). The drift time measuring algo- 
rithm was then applied to the simulated pulses so that a comparison with 
the drift time obtained from the simulation directly could be made. 

Panel (a) of Figure [3] shows the correlation between the drift time mea- 
sured using the method described above and that obtained from the simu- 
lation for a set of simulated pulses with 1% noise (relative to the amplitude 
of the pulses) distributed within the outer 5 mm of the radius of the crystal. 
As the drift time becomes longer, the measured drift time diverges as the 
initial part of the pulse becomes difficult to distinguish from noise. Panel (b) 
shows that the simulated drift time is directly proportional to the z position 
of the charge depositions in the crystal, where z = corresponds to the point 
contact. Panel (c) shows how the measured drift time correlates with posi- 
tion in the crystal. The divergence starts for events occurring in the ~35% 
of the volume furthest from the point contact. For events in the remaining 
65% of the volume closest to the point contact, the measured drift time and 
that from the simulation are very well correlated. The data simulated with 
0.1% noise shows a smaller divergence since the noise has a smaller effect. 
Including pulses from the bulk of the crystal widens these distributions only 
slightly as the drift in the z direction contributes more to the total drift time 
than the drift in the radial direction. 



4. Data collection 

Data were collected with a custom PPC HPGe detector. The 61.5 mm 
(dia.) by 50 mm crystal was purchased from ORTEC [13] and fabricated into 
a detector at Lawrence Berkeley National Laboratory (LBNL). The point 
contact is a hemispherical dimple with a diameter of 1mm ± 0.5 mm. The 
crystal has a lithium-diffused n-|- contact with aluminum metalization that 
covers the entire surface, except for the side with the point contact where 
the lithium- diffused region only wraps around for the outer 5 mm and is not 
aluminized. A schematic is shown in Figure HI The crystal was held in 
a prototype copper mount designed for the Majorana Demonstrator 



experiment [3|, ll3|, [ij] . A custom, low-radioactivity electronic board with a 
FET that can be operated at cryogenic temperature was mounted close to 
the point contact, and the amplification loop was closed in a custom charge- 
sensitive preamplifier outside the cryostat. Processed pulses from the charge 
sensitive preamplifier were read out with a 16-bit 100 MHz Struck SIS3302 
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Figure 3: Panel (a) shows the correlation between the measured drift time (y-axis) and 
aetual drift time (x-axis) for simulated pulses with 1% Gaussian noise (similar to the noise 
in 60keV events). Panel (b) shows how the actual drift time correlates with the Z-position 
of interactions in the crystal. Panel (c) shows how the measured drift time varies with 
position in the crystal. The measurement of drift time diverges from a linear relation at 
longer drift times; these events correspond to charge depositions far from the point contact 
and the full drift time become difficult to measure because of noise. About 35% of events 
distributed uniformly in volume end up in this region. 

(firmware version 1408) digitizing card. This detector and the associated 
electronics will be discussed in a later publication. 

5. Position dependence of drift time 

The drift time measurement was made using data from a coUimated Am 
source directed at different axial positions (z axis) on the crystal. The ^^^Am 
source emits 59.5 keV gamma rays that have a scattering length of order 1 mm 
in germanium so that the position of the collimated source is well correlated 
with the position of the interactions in the crystal. Figure shows the mean 
drift time of full energy events as a function of the source position. The drift 
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Figure 4: Side view schematic of PPC detector fabricated at LBNL, showing the point 
contact (A), the lithium-diffused n+ contact (B and C), and the vapor-deposited ahiminum 
ohmic contact for the high voltage (D). The thickness of the lithium-diffused and aluminum 
layers are not to scale. The z axis referred to in the text points up with the origin at the 
point contact. 

time for each pulse was determined as tgo — to.i, using the algorithm from 
equation [H with a 500 ns averaging time (50 samples for a 100 MHz digitiza- 
tion). The mean drift time was obtained for each position of the source by 
fitting the distribution of drift time of many (~ 1000) events between 58 keV 
and 62keV using a Gaussian distribution. The error on the mean drift was 
taken as the standard deviation of the fitted Gaussian. The dependence of 
the mean drift time as a function of z is consistent with that observed in the 
simulated pulses. This indicates that, for the ~ 30 mm closest to the point 
contact, the measured drift time varies linearly with the axial position in the 
crystal. 

One can thus use the distribution of drift times from a detector to mon- 
itor the homogeneity of interactions in the crystal. This can be of use for 
low background experiments, where contaminated parts could be sources of 
background. An optimized arrangement of PPC crystals could, in principle, 
be used in an array to monitor the isotropy of external backgrounds. Addi- 
tional sensing electrodes could be added to this type of detector to enhance 
the precision of the interaction sites in the crystal and lead to gamma ray 
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tracking applications. 
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Figure 5: Mean drift time of events as a function of z position of a collimated ^''^ Am source. 
For each source position, the distribution of measured drift times of events between 58keV 
and 62 keV was fit to a Gaussian. The mean and standard deviation from the fit were 
used to define the mean drift time and error as a function of z. The trend compares well 
to panel (c) of Figure |3] which was done for the simulated pulses. 



6. Charge trapping model 

Following jsj one can use the measured drift time of a cloud of charge 
carriers to examine the effects of their trapping by impurities in the crystal by 
implementing a simple model. The effective trapping length for this process 
can be modeled in the time domain as a "charge-trapping time" , r, related to 
the charge carrier lifetime in the detector. As the drift time, tdrift, increases, 
the amount of charge in the cloud, Q{tdrift) decreases exponentially from the 
initial deposited charge, Qo- 

^drift 

Qitdrift) = Qoe - (2) 

Since the energy deposited is proportional to the charge in the cloud, a loss of 
charge due to trapping results in an underestimate of the deposited energy 
if no correction is applied. Given the uncorrected measured energy, Emes, 
and the charge trapping time, the corrected measure of the deposited energy, 
Ecorr, can be recovered: 



9 



7. Charge trapping characterization and energy resolution improve- 
ment 

In order to examine charge trapping effects, ^"^Co and ^^^Th source data 
were collected using the same detector setup as described above. Of par- 
ticular interest is the double escape peak (DEP) at 1592 keV from the ^^^Tl 
2614 keV line in the ^^^Th decay chain. Events in the DEP are intrinsically 
"single-site events" (SSE); that is, over 95% of these events have only one en- 
ergy deposition in the crystal and allow one to calibrate methods of rejecting 
multi-site events. Indeed, one would expect the charge trapping correction 
described above to work best for events having only 1 interaction in the 
crystal, as multi-site interactions would contain multiple components of drift 
time. Conversely, events from the single escape peak (SEP) at 2104 keV are 
intrinsically multi-site. 

Multi-site events can effectively be rejected in PPC detectors using the 
"A/E method' 0, [6|. This method compares the maximum amplitude of the 
current pulse, A (calculated using equation [1] with a shorter averaging time 
of 200 ns), to the calibrated energy, of the event. In a multi-site event, the 
current pulse will have multiple peaks, whereas in a single-site interaction, 
the maximum amplitude of the current pulse will be proportional to the 
deposited energy. Charge trapping effects are examined with and without the 
A/E-based selection of single-site events. The crystal used in this experiment 
was shown to reject 93% of multi-site events while preserving 90% of single- 
site events using the A/E method. 

Charge trapping effects in a detector can be characterized by applying 
equation [3] to the measured energy of events and determining the energy 
resolution as the charge trapping time, r, is varied. The drift time defined 
in section [2] is used to apply equation [3] for each event. Figure [6] shows the 
full width at half maximum (FWHM) energy resolution of the 1332 keV ^°Co 
line as r is varied, with and without the SSE-selection cut applied. Plotting 
these curves is a good tool for examining charge trapping effects in a crystal. 
If no charge trapping occurred, there would not be a global minimum since 
increasing the charge trapping time indefinitely is equivalent to applying no 
correction. The FWHM resolution with no correction is 2.30 keV for the raw 
spectrum and 2.60 keV for the SSE-selected spectrum; the asymptotic value 
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is reached beyond the range plotted in the figure. Using these curves, the 
optimal charge trapping time is determined by fitting a parabola to the data 
in the region of the minimum. The uncertainty on r is determined using the 
uncertainty on the energy resolution at the minimum of the fitted parabola; 
the range in r where the parabola is within the uncertainty on the resolution 
is taken as the error on the determined value of r. 

It is interesting to note that, while the application of the single-site event 
selection cut degrades the energy resolution (from 2.30 keV to 2.60 keV, be- 
yond the range shown in Figure |6]), the charge trapping correction removes 
this effect. This can be explained by the fact that after the SSE events 
are selected, the sample of events contains a larger fraction of events with 
large drift times. The multi-site events with large measured drift times con- 
tain a significant contribution from the difference in the arrival times from 
multiple charge depositions rather than from actually long drift times. For 
this reason, a sample of events selected as single-site is more susceptible to 
charge trapping. This is an important effect to understand if one applies 
a single-site selection cut to remove backgrounds (such as in a neutrinoless 
double-beta decay experiment) and were to naively extrapolate the energy 
resolution that is measured using (predominantly) multiple-site events from 
calibration data. 

Figure [7] shows the resolution of the 1332 keV peak as different start times 
are used to define the drift time. There is a substantial improvement in energy 
resolution when using earlier start times to define the drift time. Incidentally, 
the start time defined using the 10% point of the average current is similar 
to that obtained from using the 10% rise of the charge pulse, as used in jsj. 
This figure is indicative of the improvement that is obtained from using an 
optimized definition of drift time. 

Using peaks at 727keV, 911keV, 1173 keV, 1332 keV, 1592 keV and 2104 keV, 
the optimal charge trapping time and uncertainty were found for each energy 
(and tabulated in tables [1] and [2]). An average charge trapping time was 
determined by weighing the value at each energy with its uncertainty and was 
then used to correct the entire energy spectrum. The average charge trapping 
time for this crystal was determined to be ~ 460^s. Figured shows the energy 
spectrum before (solid) and after (dashed) the charge trapping correction was 
applied, for events with and without the SSE-selections. Panels (a) and (b) 
show the energy spectrum over a large range, whereas panels (c) and (d) show 
the spectrum in the region around the "^^^Tl DEP, highlighting the effect of 
the SSE-selection. In particular, one notes the substantial reduction in the 
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Figure 6: Energy resolution of the ^°Co 1332 keV peak as a function of charge trapping 
time with SSE-selection applied (horizontal hash, light line) and for the raw spectrum 
(vertical hash, dark line) . The region near the photopeak was fitted to a Gaussian with a 
linear background. The standard deviation of the fit and its associated fit error were used 
to determine the FWHM and the resulting band. The error band for the data with SSE- 
selection applied is larger because of statistics (since most 1332 keV events are multi-site 
and are thus removed). The presence of a global minimum (at r w 450/is) is indicative 
of charge trapping in the crystal. As the charge trapping time increases, the resolution 
asymptotically tends to its uncorrected value. 

^^^Ac peak at 1588 keV. The charge trapping correction introduces a change 
in energy scale that was accounted for in these figures; it was verified that 
this change in energy scale was linear (see equation [3]) . 

The improvement in energy resolution from the charge trapping correc- 
tion is summarized in Tables [1] and [2] for all events and for those selected 
as single-site, respectively. The results presented here show that the charge 
trapping correction originally suggested by jsj greatly benefits from an op- 
timal determination of the drift time. For all energies, the charge trapping 
correction removes the bias in energy resolution introduced from applying the 
SSE-selection; the corrected energy resolution is now consistent between all 
events and those selected as single-site. The data presented here show that 
charge trapping corrections can be important for ^^Ge neutrinoless double- 
beta decay experiments that strive to have good energy resolution while 
implementing cuts to remove multiple-site interactions. 
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Figure 7: Energy resolution of the 1332 keV peak as a function of the start time that was 
used to determined the beginning of the drift time, as measured in percent of the calculated 
average current waveform. The charge trapping correction is improved dramatically when 
using earlier start times. 
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Table 1: Summary of charge trapping correction performance for selected lines from Figure 
UK (no SSE-selection applied). The table shows the nominal FWHM energy resolution 
(with no correction), the best resolution found by varying t, the value of r that gives the 
best resolution, the resolution at the average value of r and the relative improvement in 
energy resolution at the average value of r for various lines in the energy spectrum for 
^''Co and isotopes in the ^■^^Th decay chain. Uncertainties are shown in parentheses. 

8. Conclusion 

With PPC detectors gaining popularity as the detector of choice in various 
types of low background experiments using germanium, understanding their 
properties has become an area of increased interest. This work presented a 
method to estimate the drift time of charge carriers in PPC detectors. Given 
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Table 2: Summary of charge trapping correction performance for selected lines from Figure 
[5}d (events selected as single-site). The table shows the nominal FWHM energy resolution 
(with no correction), the best resolution found by varying t, the value of r that gives the 
best resolution, the resolution at the average value of r and the relative improvement in 
energy resolution at the average value of r for various lines in the energy spectrum for 
®°Co and isotopes in the ^'^^Th decay chain. Uncertainties are shown in parentheses. The 
^°^T1 single escape peak is not included since the SSE cut removes essentially all of those 
multi-site events. 

the drift time, one has some sensitivity in determining the position of an 
energy deposition, which can be of use for monitoring the isotropy of events 
in a low background detector. It has also been shown, following previous 
work by [sl, that one can use the drift time to look for evidence of charge 
trapping in PPC detectors. A correction to improve the energy resolution of 
PPC detectors was found to benefit greatly from using an optimized mea- 
sure of drift time. It was noted that selecting single-site events using the 
A/E method can result in an apparent increase in charge trapping effects 
and a biased estimate of the energy resolution. Since the energy resolution 
has a direct impact on the sensitivity of neutrinoless double-beta decay ex- 
perimental searches, these types of experiments must correctly account for 
charge trapping effects when determining their energy resolution if they use 
an A/E- type selection to reduce backgrounds. 
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Figure 8: Energy spectrum of events before (solid) and after (dashed) a eorreetion for 
charge trapping was apphed. Panels (a) and (b) show a larger range of energies, whereas 
panels (c) and (d) highlight the region around the ^^^Tl DEP at 1592 keV. Spectra are 
shown for all events (panels (a) and (c)) as well when the SSE-selection is applied (panels 
(b) and (d)). A change in energy scale from the charge trapping correction was accounted 
for in these figures. 
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